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in the plane perpendicular to the magnetic field. The relevance of these results for cosmological 
models based on supersymmetric dynamics is pointed out. 
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I. INTRODUCTION. 



> _ 

It has been recently shown that a constant magnetic field in 2 + 1 and 3 + 1 dimensions is a strong catalyst 
of dynamical chiral symmetry breaking, leading to the generation of a fermion dynamical mass even at the weakest 
I attractive interaction between fermions. 

The essence of the effect is the dimensional reduction D ^ D — 2 (i.e. 2+1^0+1 and 3 + 1 — > 1 + 1 ) 
in the infrared dynamics of the fermion pairing in a magnetic field. The physical reason of this reduction is the 
Q\ , fact that the motion of charged particles is restricted in those directions that are perpendicular to the magnetic 
field. This is in turn connected with the point that, at weak coupling between fermions, the fermion pairing, leading 
to the chiral condensate, is mostly provided by fermions from the lowest Landau level (LLL) whose dynamics are 
ft, ' {D — 2)-dimensional. 

I In this paper, we shall further clarify the effect of the dynamical reduction, studying in detail the infrared dynamics 
^*| 1 in the (3 + l)-dimensional Nambu-Jona-Lasinio (NJL) model in a magnetic field. We shall consider both the ordinary 
(non-supersymmetric) and supersymmetric versions of the NJL model. In particular, we shall show that in the 
• i-H "continuum" limit, when both the strength of the magnetic field and the ultraviolet cutoff go to infinity, both the 
non-supersymmetric and supersymmetric (weakly coupling) NJL models with iV c colors are reduced to a continuum set 
of independent (l+l)-dimensional Gross-Neveu (GN) models ||, labeled by coordinates xj_ in the plane perpendicular 
to the magnetic field B. The number of colors in the GN models is N c = (ir/2C)N c , where C is C = A 2 /\eB\ in the 
"continuum" limit (here A is the ultraviolet cutoff). As will be shown in Sec. 3, the factor tt/2C is proportional to a 
(local) magnetic flux attached to each point in the a;j_-plane. 

On the other hand, at strong coupling, the dynamics in the supersymmetric and non-supersymmtric NJL models 
are very different. 

Recall that there is no spontaneous chiral symmetry breaking in the supersymmetric NJL model ||. An external 
magnetic field changes the situation dramatically: chiral symmetry breaking occurs for any value of the coupling 
constant in this model. This agrees with the general conclusion of Refs. J|-§| that the effect of the catalysis of chiral 
symmetry breaking by a magnetic field in 3 + 1 dimensions is a universal, model-independent effect. 

As was already shown in Ref. , the dimensional reduction 3 + 1— > 1 + 1 in the dynamics of the fermion pairing in 
a (finite) magnetic field is consistent with spontaneous symmetry breaking. Recall that, due to the Mermin-Wagner- 
Coleman (MWC) theorem 0, there cannot be spontaneous breakdown of continuous symmetries at D = 1 + 1. The 
MWC theorem is based on the fact that gapless Nambu-Goldstone (NG) bosons cannot exist in 1 + 1 dimensions. 
However, in a magnetic field, the reduction 3 + 1 — > 1 + 1 takes place (in the infrared region) only for propagators of 
charged particles: it reflects the fact that the motion of charged particles is restricted in the directions perpendicular 
to the magnetic field. On the other hand, NG bosons connected with spontaneous chiral symmetry breaking are 
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neutral and therefore their propagators have (3 + l)-dimensional form QR. This in turn implies that the effect of 
spontaneous chiral symmetry breaking does not contradict the MWC theorem. 

The interplay between the GN model and the NJL model in a magnetic field established in this paper further 
clarifies this issue. As was mentioned above, the NJL model is reduced to a set of the GN models only as the strength 
of the magnetic field goes to infinity. As we shall show in Sec. 4, at finite \eB\, the dynamics of the NJL model in 
a magnetic field are in a sense similar to the dynamics in the (2 + e)-dimensional GN model: the magnetic length 
I = |ei?| _1 / 2 plays here the role of the (physical) e-parameter which is an infrared regulator. 

As was already pointed out in Refs. |l]-f|], there may be interesting applications of this effect in cosmology as 
well as in particle and condensed matter physics. The results of the present paper may be particularly relevant for 
cosmological scenarios based on supersymmetric dynamics |J . 

The paper is organized as follows. In Section 2 we, for completeness, derive the effective action in the GN model. 
In Sections 3 and 4 we establish the connection between the GN model and the NJL model in a magnetic field as 
\eB\ — > co. In Section 5 we consider the dynamics of the supersymmetric NJL model in a magnetic field. In Section 
6 we summarize the main results of the paper and discuss possible applications of these results and as well as the 
possibility of their extension to inhomogeneous magnetic field configurations. In the Appendix some useful formulas 
and relations are derived. 



II. EFFECTIVE ACTION IN THE GROSS-NEVEU MODEL. 



In this section, for completeness, we shall derive the effective action for the GN model. The Lagrangian density of 
the GN model is: 

Cgn = \ [*, (iyfy)*] + y [(**) 2 + (^7 5 *) 2 ] (1) 

where fj, = 0, 1 and the fermion field carries an additional "color" index a = 1, 2, . . . , N c (for simplicity, we consider 
the case of the chiral Ul(1) X Ur(1) symmetry). The theory is equivalent to the theory with the Lagrangian density 

z'gn = \ (<y^)*] - + *t 5 o* -^( (j2 + ■ ( 2 ) 

The Euler-Lagrange equations for the auxiliary fields a and tt take the form of constraints: 

a = -Gftf, tt = -Gfry 5 *, (3) 

and the Lagrangian density (|[) reproduces Eq.(Q) upon application of the constraints (^|). The effective action for the 
composite fields a and tt can be obtained by integrating over fermions in the path integral. It is given by the standard 
relation: 

r G Jv(<x, tt) = f GN (a, tt) - X, J d 2 x{a 2 + tt 2 ), (4) 

f GN {a, tt) = -iTrLn [i-y"d„ - (a + i 7 5 7r)] . (5) 

The low energy quantum dynamics are described by the path integral (with the integrand exp(iLcAr)) over the 
composite fields a and tt. As N c — > co, the path integral is dominated by the stationary points of the action: 
STgn/Sct = STgn/^ — 0. We will analyze the dynamics by using the expansion of the action Tgn in powers of 
derivatives of the composite fields. 

We begin the calculation of Ton by calculating the effective potential Vgn- Since Vgn depends only on the 
Ul(1) x [/fl(l)-invariant p 2 = a 2 + 7r 2 , it is sufficient to consider a configuration with tt = and a independent of x. 
Then we find from Eqs. (|) and (§): 



x The Lorentz invariance is broken by a magnetic field in this problem. By the (3 + l)-dimensional form, we understand 
that the denominator of the propagators depends on energy and all the components of the center-of-mass momentum, i.e. 
V{P) ~ (P 2 - C||P,? - CiPi)" 1 with C«,C± / 0. 
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where the integration is done in Euclidean region (A is an ultraviolet cutoff). As is known, in the GN model, the 
equation of motion dVcpj/dp = has a nontrivial solution p = a = md yn for any value of the coupling constant G. 
Then the potential Vgn can be rewritten as 
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Due to the MWC theorem (t]] , there cannot be spontaneous breakdown of continuous symmetries at D = 1 + 1. The 
parameter md yn is an order parameter of chiral symmetry breaking only in leading order in 1/ N c (this reflects the 
point that the MWC theorem is not applicable to systems with N c — > oo [||). In the exact GN solution, spontaneous 
chiral symmetry breaking is washed out by interactions (strong fluctuations) of would-be NG bosons tt (i.e. after 
integration over 7r and a in the path integral). The exact solution in this model presumably corresponds to the 
realization of the Berezinsky-Kosterlitz-Thouless (BKT) phase: though chiral symmetry is not broken in this phase, 
the parameter mdyn still defines the fermion mass, and the would-be NG boson 7r transforms into a BKT gapless 
excitation ||. 

Let us now turn to calculating the kinetic term in Ton- The chiral x Ur(1) symmetry implies that the 

general form of the kinetic term is 



: gn = -^-( d ^Pj d ^Pj) + -jr(Pj d nPj)(pidvPi) 
z p 



(9) 



where p = (a, it) and f^", j^" are functions of p 2 . To find the functions f^ v and f% v , one can use different methods. 
We utilize the same method as in Ref. jlj (see Appendix A in that paper). The result is: 
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with S(k) = «(fc M 7 Al + p)/(k 2 — p 2 )- The explicit form of these functions is: 
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III. THE INTERPLAY BETWEEN THE GN MODEL AND THE NJL MODEL IN A MAGNETIC FIELD 

In this section, we compare the effective actions in the GN model and in the NJL model in a magnetic field, and 
we establish a rather interesting connection between these two models. 

The analog of the Lagrangian density (||) in the NJL model in a magnetic field is 

C = \ H^D^] - *(a + i 7 5 7r)* - ± (a 2 + ir 2 ) (13) 

where D a = — ieAff*, A e * 1 — Bx 2 S^ (the magnetic field is in +Xi direction). 

In leading order in 1/N C , the effective action in the NJL model in a magnetic field is derived in Refs. [^Q. The 
effective potential and the kinetic term are (p = (a, 7r)): 
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with /" and ft, being diagonal tensors: 
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Here G is the NJL coupling constant, N c is the number of colors, Q(y,x) is the generalized Riemann zeta function, 
Cf{v,x) — d£(v,x)/dv, 7 w 0.577 is the Euler constant, VK 3 -) = d(lnT(x)) /dx, and I = |ei?| _1 / 2 is the magnetic 
length. The gap equation dV/dp = isQ 
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where the dimensionless coupling constant g = N C GA 2 /Att 2 . In the derivation of this equation, we used the relations 



— £(u,x) = -vC{v + l,x), 
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lnT(.x) - -ln27r, C(0,Jc) 



(18) 
(19) 



As B — > (I — > oo), we recover the known gap equation in the NJL model (for a review see Ref. pJ|): 
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This equation admits a nontrivial solution only if g is supercritical, g > g c — 1 (as Eq. (|l3|) implies, a solution to the 
gap equation, p = a, coincides with the fermion dynamical mass, a = rridyn)- As was shown in Refs. [|,|, at B ^ 0, 
a non-trivial solution exists for all g > 0. 

Let us consider the case of small subcritical g, g <C g c = 1, in detail. A solution is seen to exist for this case if pi is 
small. Specifically, for g <C 1, the left-hand side of Eq.(^) is positive. Since the first term of the right-hand side in 
this equation is negative, we conclude that a non-trivial solution to this equation may exist only for 



p 2 ln(A0 2 «iln-i F 



(21) 



2 In this paper we consider the case of a large ultraviolet cutoff: A 2 S> er 2 , \eB\, where a is a minimum of the potential V. 
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(T(p 2 l 2 /2) « 2/{pl) 2 ). We then find the solution: 
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Actually, since Eq.(p^) implies that condition ( pi] ) is violated only if (1 — g) < |ei?|/A 2 , the expression 
for all 5 outside that (scaling) region near the critical value g c = 1. Note that in the scaling region (<? — 
expression for m d is different 0,0]: 
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Let us compare relation (g2|) with relation (g) for the dynamical mass in the GN model. The similarity between them 
is evident: \eB\ and |e£?|G in Eq.(p2|) play the role of an ultraviolet cutoff and the dimensionless coupling constant 
G in Eq.(^J). Let us discuss this connection and show that it is intimately connected with the dimensional reduction 
3 + 1— ► 1 + 1 in the dynamics of the fermion pairing in a magnetic field. 

Eq.(|[) implies that the GN model is asymptotically free, with the bare coupling constant G — 2ir/N c ln(A 2 /mjiL n ) — > 
as A — > co. Let us now consider the following limit in the NJL model in a magnetic field: \eB\ — > oo, A 2 /\eB\ = 
C 1. Then relation (E2), which can be rewritten as 
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implies that the behavior of the bare coupling constant g must be 
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in order to get a finite value for ra 2 dyn in this limit. 



Thus in this "continuum" limit, we recover the same behavior for 



the coupling g in the NJL model as for the coupling constant G in the GN model. 

Let us now compare the effective potentials in these two models. At first glance, the expressions (§ and (0) for 
the effective potentials in these models look very different: the character of ultraviolet divergences in 1+1 and 3+1 
dimensional theories is essentially different. However, using Eqs.(|l8|) and ([l9]), the expression ( |l4| ) can be rewritten, 
for small pi, as 
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Then, expressing the coupling constant g through from Eq.(|2^), we find that 
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Here we used the fact that, because of Eq. (pi|) , the ratio (pi) 2 is small near the minimum p = m c i yn - 

The expressions (|7) and (27) for the potentials in these two models look now quite similar. There is however an 

additional factor \eB /2tt in the expression (p7j). Moreover, the field p, which depends on the two coordinates Xq and 

x\ in the GN model, depends on the four coordinates xo, x\, x-i and X3 in the NJL model. 

In order to clarify this point, let us turn to the analysis of the kinetic term ([l5]) in the effective action of the NJL 

model in a magnetic field. 

Because of the expression (E2) for mdyn at small g, the term l/(pl) 2 dominates in the functions /j 30 = — ff 1 and 
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Up to the additional factor \eB\/2-K, these functions coincide with those in (12) in the GN model. On the other hand, 
the functions / 22 = /f 3 and /f 2 = /| 3 , connected with derivatives with respect to the transverse coordinates, are 



strongly supressed, as compared to the functions (|28|), and the ratios of the functions f± = fi and f% = f% to 
those in (g8|) go rapidly (as rn 2 dyr J A 2 ) to zero as \eB\ — > oo, A 2 /\eB\ = C. 

As a result, the coordinates x 2 and 23 become redundant variables in this limit: there are no transitions of 
field quanta between different points in the £2X3— plane. Therefore the model degenerates into a set of independent 
(1 + l)~dimensional models, labeled by xi = (x%, X3) coordinates. Let us show that all these models coincide with 
a (1 + 1) -dimensional GN model with the number of colors N c = (tt/2C)N c , where the factor C = A 2 /\eB\ was 
introduced above, in the definition of the "continuum" limit (\eB\ — > 00, A 2 /\eB\ — C). 

Let us put the N JL model on a lattice with a the lattice spacing of the discretized space-time (in Euclidean region) . 
Then its effective action can be written as 

F N jl(ct, 7r) = / dx 2 dx 3 / dx i dxiL { -^/^{a(x),iT{x)) 



1 ^ ^ 

-^ eB ^ E E ^JLK(n»"»),7r ij -(n,m)) (29) 

i 7 j= — OO 71,771= — OO 

where a^ (n,m) — <r(x), 7r^(n,m) = 7r(x), with X2 — ia, ^3 — ja, %a = i%o = na, X\ = ma, and here the factor 
\eB\/2it was explicitly factorized from L^/l • Now, taking into account Eqs.(0),(|l^) and Eqs.(^), (^9|), we find that 

Tnjl = 2^' eB ' a4 E E ^vjiVyK m ),^jK™)) 

i,j= — 00 n,m— — 00 

-> E / dx i dx ^ L GN ] ( a x 2 x 3 (X\\), Kx 2 x 3 (X||)) (30) 

in the "continuum" limit with (|e£?|/27r)a 2 = 7r|eS|/2A 2 = ir/2C (here A = 7r/a is the ultraviolet cutoff on the 
lattice)^. The lagrangian density corresponds to the GN model with the number of colors iV c = (ir/2C)N c . 

Note that the symbol ^ here is somewhat formal and it just implies that the GN model occurs at each point in 

the X2X3-plane. 

The physical meaning of this reduction of the NJL model in a magnetic field is rather clear. At weak coupling, the 
fermion pairing in a magnetic field takes place essentially for fermions in the LLL with the momentum k\=0. The size 
of the radius of the LLL orbit is l=\eB\~ x l 2 [[l2[ . As the magnetic field goes to infinity, this radius shrinks to zero. 
Then, because of the degeneracy in the LLL ]12|, there are (|e£?|/27r)a 2 = n/2C states with ki—0 at each point in the 
X2X3~plane. This degeneracy factor (proportional to the (local) magnetic flux across a plaquette) leads to changing 
the number of colors, N c — + iV c = (ir/2C)N c , in the GN model. Note that since N c appears analytically in the path 
integral of the theory, one can give a non-perturbative meaning to the theory with non- integral N c . 

A few comments are in order. 

Since these GN models are independent, the parameters of the chiral Ul(1) x Ur(1) transformations can depend 
on x_i_. In other words, here the chiral group is J! U { L x±) (l) x U£ ±} (1). As a result, there are an infinite number of 

x± 

gapless modes ^^^(xn) in the "continuum" limit. 

Since there is no spontaneous breakdown of continuous symmetries at D = 1 + 1, the fields Tr^^^ll) do not describe 
NG bosons (though they do describe gapless BKT excitations) || . 

Since the magnetic field depends on A in the "continuum" limit, it can be considered as an additional parameter 
("coupling constant") in the renormalization group. The ratio b = \eB\/A 2 = C _1 is arbitrary here. ^From the point 
of view of the renormalization group, this can be interpreted as the presence of a line of ultraviolet fixed points for 
the dimensionless coupling b. The values of b on the line define the local magnetic flux and, therefore, the number of 
colors N c in the corresponding GN models. 



3 Of course, the ultraviolet cutoff on the lattice is different from the cutoff in the proper-time regularization used above. 
However, since the constant C = A 2 /\eB\ is anyway arbitrary here, we use the same notation for the cutoff on the lattice as in 
the proper-time regularization. 



G 



Our consideration of reducing the NJL model in a magnetic field to a continuum set of the GN models was somewhat 
heuristic. It would be worth deriving this reduction in a more rigorous way, putting the NJL model on a lattice in 
Euclidean space and then realizing renormalizations in the "continuum" limit. 

In the next section, we shall discuss the connection between the GN model and the NJL model in a magnetic field 
in more detail. 



IV. MORE ABOUT THE CONNECTION BETWEEN THE GN MODEL AND THE NJL MODEL IN A 

MAGNETIC FIELD 



In the previous section, the reduction of the NJL model in a magnetic field of the infinite strength to a continuum 
set of the GN models was established. But what is the connection between the NJL and GN models at finite, though 
large, values of the magnetic field? 

In order to answer this question, let us turn to a more detailed discussion of the infrared dynamics within these 
two models. 

The GN model is asymptotically free, with the bare coupling constant G — 2ir/N c \n(A 2 /m 2 yn ) — > as A — > oo. 
Therefore there is dimensional transmutation in the model: in the scaling region (G <§C 1) the infrared dynamics, 
with momenta k satisfying [ln(A 2 //c 2 )] _1 -C 1, are essentially independent of either the coupling constant G or the 
cutoff A; the only relevant parameter is the dynamical mass m^yn (which is an analogue of the parameter Aqcd m 
QCD). Because of Eq.(22) for mdyn, one might expect that a similar dimensional transmutation should take place in 
the NJL model in a magnetic field: for \eB\GN c <C 1, the infrared dynamics, with k 2 <§; \eB\, should be essentially 
independent of the magnetic field and the coupling G. 

The real situation, however, is more subtle. Let us look at the propagators for fermions and a and ir particles at 
low momenta in the NJL model with a magnetic field, in leading order in 1/N C M: 
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(31) 
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where = (fcg, ^3) and where C„ and C w are some inessential constants (we consider the case of the weak coupling, 
when relation (|22| ) is valid). Because of relations (22) and (p4|), the coefficients of the k^-terms in the propagators 
T> a {k) and T>^(k) are exponentially small. Also, in the infrared region, the fermion propagator is independent (up 
to power corrections ~ (kj_/|e.B|) n , n > 1) of the magnetic field. Therefore one might think that in this case, like 
in the GN model, the only relevant parameter for the infrared dynamics is mdyn- However, while the dependence of 
the propagators of fermions and a particles on \eB\ can indeed be neglected, this dependence is essential in the case 
of the propagator of the gapless NG mode ir. The k^-term in T> n (k) provides the (3 + l)-dimensional character of 
the NG propagator which, as explained in Introduction, is crucial for the realization of spontaneous chiral symmetry 
breaking in the model. In a sense, the magnetic length I = \eB\~ x / 2 plays here the same role as the e-parameter in 
the (2 + e)-dimensional GN model ||. 

As \eB\ — > 00, the NJL model is reduced to a set of independent GN models. However, at finite values of \eB\, 



the transverse velocity of the NG mode, though small [|vj_| < y (TO 2 yn /|eB|) ln(|eS|/7rrn 2 yn )], is not zero. Therefore 

there are now transitions of field quanta between different points in a^a^-plane, i.e. interactions occur between the 
GN models associated with different points xj_. As a result, at finite \eB\ the chiral group is reduced to global 
C/l(1) x Ur(1) transformations, and there is only one NG boson tt for this case. 

Therefore the reduction of the NJL model, described in the previous section, takes place only as \eB\ ^00. At 
finite values of the magnetic field, the dynamics in the NJL and GN models are different: while there is spontaneous 
chiral symmetry breaking in the NJL model, the BKT phase is realized in the GN model [|| . The connection between 
these two sets of dynamics is similar to that between the dynamics of 2-dimensional and (2 + e)-dimensional GN 
models. 

In conclusion, we emphasize that this discussion pertains only to the NJL model with a weak coupling constant, 
when relation ( |22j ) is valid. In the case of the NJL model with a near-critical g, the situation is different: when 
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g — > g c — 0, ( p3| ) is valid. The difference between these two dynamical regimes reflects the fact that, while at weak 
coupling the LLL dominates, at near-critical g all Landau levels are relevant [Q . 

In the next section, we shall discuss the dynamics in the supersymmetric NJL model. 

V. CHIRAL SYMMETRY BREAKING IN THE SUPERSYMMETRIC NJL MODEL IN A MAGNETIC 

FIELD 

As is well known, there is no spontaneous chiral symmetry breaking in the supersymmetric (SUSY) NJL model 
at any value of the coupling constant |^|. Here we shall show that an external magnetic field changes the situation 
dramatically: chiral symmetry breaking occurs for all (positive) values of the coupling constant. Moreover, at weak 
coupling, the dynamics in SUSY and ordinary NJL models are similar (and, therefore, intimately connected with the 
dynamics in the (1 + I)-dimensional GN model). 

The action of the SUSY NJL model with the Ul(1) x Ur(1) chiral symmetry in a magnetic field is 



r = J d 8 z[Qe v Q + Q c e- v Q c + G{Q C Q){QQ C )] 



(34) 



Here we utilize the notations of Ref. [[L3|, except for our choice of metric g^ v — diag(l, — 1, — 1, — 1). In Eq. (p4[) , 
d 8 z = d 4 xd 2 9d 2 9, Q a and Q c a are chiral superfields carrying the color index a — 1,2,... ,N C , i.e. Q a and Q c a are 
assigned to the fundamental and antifundamental representations of the SU(N C ), respectively: 



Q a = <p a + V29^ a + 9 2 F a , Q c a = <p c a + V29i/> c a + 9 2 F^ 



(35) 



(henceforth we will omit color indices). The vector superfield V{x,9,9) = —9a tl 9A e ^ t , with A?, xt = Bx 2 S^, describes 
an external magnetic field which is in the +x\ direction. 
The action (134) is equivalent to the following action: 



r A = / d 8 z 



Qe v Q + Q c e~ v Q c + ^HH 



G 



d (i : 



^HS - QQ C S 



d 6 z 



-HS-QQ C S 
G 



Here d 6 z = d 4 xd 2 9, d 6 z = d 4 xd 2 9, and H and S are two auxiliary chiral fields: 

H = h + V29 Xh + 9 2 f h , S = s + V29 Xs + 9 2 f s . 
The Euler-Lagrange equations for these auxiliary fields take the form of constraints: 



H = GQQ C , S = --D 2 (H) 



G 



D 2 {QQ C ). 



(36) 



(37) 



(38) 



Here D is a SUSY covariant derivative |13j. The action ( pq ) reproduces Eq.(34) upon application of the constraints 
In terms of the component fields, the action (|36|) is 



d^x 



+ F C 'F C + - {-tfOh + iXh^d^Xh + flfh 
+ ^ (XhXs - hf 8 - sf h + h.c.) 

- (Wv c + (^ c + f c ^)Xs - s{ V F c + if c F) - W C U + h.c) 



(39) 



Let us consider the effective potential in this model. For this purpose, one can treat all the auxiliary scalar fields as 
(independent of x) constants and all the auxiliary fermion fields equal zero (since the auxiliary fields are colorless, 
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loop diagrams involving them do not contribute to the effective potential in leading order in X/N c ). Then, using the 
Euler-Lagrange equations for the fields F, F c , fh, h and their conjugates, we find = —s(p c , = —sip, fl = s, 
fj = 0, plus h.c. equations. Then the action becomes 



d 4 x 



J [(d, - ieA^f +p 2 ]<p- ^ [(d, + ieAff + p 2 } v 



(40) 



where s = a + in, p 2 = \s\ 2 = a 2 + ir 2 , and the Dirac fermion field ipu is introduced. 

In leading order in 1/N C , the effective potential V{p) can now be derived in the same way as in the ordinary NJL 
model (J]. The difference is that, besides fermions, the two scalar fields ip c and <p give a contribution to V(p): 



V(p) = ^ + V fer (p) + 2V bos (p), 



where 



Vfer(p) 



8tt 2 / 4 



oo 

J -j exp (— s{lp) 2 ) coth s, 



Vbos{p) = - 



V('A) 2 

OO 

N c r ds 

16tt 2 Z 4 J ^ 

l/(/A) 



— exp i 



-s{lp) 2 ) 



I 



sinh s 



(41) 

(42) 
(43) 



As is shown in the Appendix, the potential V(p) can be rewritten as 

N, 



V{ P ) 



8tt 2 Z 4 

N c 
16tt 2 Z 4 



M + (w ^ 1 _ 1i ,<w 2 



fJ 



[('p) 2 + l]/2 

4- J dx\nT(x) 

(lpf/2 



+ 



ln(A0 2 - 7 - ln(87r 2 ) 



"'1 



(44) 



where the dimensionless coupling constant is g = GN c /8tt 2 1 2 . 

As the magnetic field B goes to zero (I — > oo), we recover the known expression for the potential in the SUSY NJL 
model Hi: 



V(p) 



El 

G ' 



(45) 



The potential V(p) is positive-definite, as has to be in a supersymmetric theory. The only minimum of this potential 
is p — corresponding to the chiral symmetric vacuum. 

The presence of a magnetic field changes this situation dramatically: at B ^ 0, a non-trivial global minimum, 
corresponding to spontaneous chiral symmetry breaking, exists for all g > 0. Moreover, we show below that, as the 
coupling g — > 0, the SUSY NJL model becomes equivalent to the ordinary NJL model and, therefore, at weak coupling 
and as B — > oo, it is reduced to the same continuum set of the (1 + l)-dimensional GN models. On the other hand, 
the dynamics of the SUSY and non-SUSY NJL models in a magnetic field are very different at strong coupling. 

The gap equation dV/dp = 0, following from Eq.(H), is 



Np 

4^2 



In 



21nT 



1 



21nL 



aw 



\ 2 



(46) 



As can be seen from (|44|), at B ^ the trivial solution p = to this equation corresponds to a maximum of V: 
d 2 V/dp 2 \ p= o = — oo. Numerical analysis of equation ( fig ) for g > and B ^ shows that there is a nontrivial solution 
p = a — rridyn which is the global minimum of the potential. The analytic expression for m^ yn can be obtained at 
small g (when m^yn I < 1) and very large g (when Tfidyn 

l » 1). In those two cases, the results are: 
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a) g <C 1 (rridynl <!)■ The gap equation (^) is 



In 



(47) 



2\eB\ 



exp 



1 

2ff 



2|eS| 



■ exp 



b) g~^> I (mdynl ^> 1)- Now, the gap equation (f46| ) gives 

1 1 



(48) 



(49) 



g\eB\ \eB\ 



47T 



(50) 



The numerical solution to the gap equation (|46|) for general values of g is shown in Fig. 1. 

Let us discuss the case of the weak coupling in more detail. By using the asymptotic series for T(x) the 
potential ( f44| ) can be rewritten at pi <C 1 as 



V(p) = V(0) 



8tt 2 Z 4 



+ (piy In 



(Plf 



- 1 



-(pl) 2 lnn + 0((pl) 4 ) 



(51) 



Using Eq.(p7|), the coupling constant g can be expressed through md yn , and the potential can be rewritten as 



V(p) = 7(0) 



8tt 2 / 2 



In- 



dyn 



i + O((p0 2 ) 



(52) 



This potential coincides (up to exponentially small terms) with that in the ordinary (weakly coupling) NJL model 
(see Eq.©). 

Now, let us consider the kinetic term in the SUSY NJL model in a magnetic field. As shown in the Appendix, at 
weak coupling, this term also coincides with that in the ordinary NJL model. Therefore, at weak coupling, these two 
models are equivalent. 

Note that the equivalence between the SUSY and non-SUSY NJL models becomes explicit only after the renor- 
malization of the coupling constants. The structure of the ultraviolet divergences in the effective potentials of these 
models is quite different (compare Eqs. Jl4"| ) and p4])). This reflects the point that only the infrared (and not ul- 
traviolet) dynamics in these models are equivalent, i.e., in the renormalization group language, these two models are 
assigned to the same universality class [[14j. 

The physical picture underlying this equivalence is clear. An external magnetic field explicitly breaks supersym- 
metry: the spectra of charged free fermions and bosons in a magnetic field are essentially different |f[^| . While for 
fermions the spectrum is 



2\eB\n + kf, 



0,1,2,. 



for bosons it is 



£„(fci) = ±Jm 2 + \eB\(2n + 1) + kf, 



0,1,2,. 



(53) 



(54) 



The crucial difference between them is that while for fermions the energy of the LLL is independent of B and (at 
k\ = 0) goes to zero as m — > 0, for bosons the energy of the LLL is E'o(fci) = ±-\/|ei?| + k\ . In other words, there is 
a gap AE = \/\eB\ in the spectrum of massless bosons in a magnetic field. Recall that, at weak coupling, the chiral 
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condensate in the NJL model occurs because of the fermion pairing in the LLL with k\ = 0. This also happens in the 
SUSY NJL model, in which the bosonic degrees of freedom become irrelevant at weak coupling. 

Notice that, in the effective action, the field p plays the role of the mass. In particular, the infrared singularities 
(as p — > 0) in the expressions for the potential and the kinetic term (see Eqs.([27|) and (28)) reflect the absence of a 
gap in the LLL for massless fermions. 

We emphasize that the equivalence between these two models takes place at weak coupling only. When g becomes 
larger than g c = 1, the dynamics in these models are essentially different. While at g > g c = 1 in the ordinary NJL 
model, the chiral symmetry is spontaneously broken even without an external magnetic field, there is no spontaneous 
chiral symmetry breaking in the SUSY NJL model as B — > at any value of g (see Eq. (|50|) ) . The difference between 
the dynamical regimes with weak and strong coupling corresponds to the fact that, while at weak g the dynamics of 
the LLL of fermions dominates, at g > g c = 1 all (fermionic and bosonic) Landau levels become relevant. 



VI. CONCLUSION. 



In this paper we have studied the infrared dynamics of both ordinary and SUSY NJL models in a magnetic field. It 
has been shown that, at weak coupling, the infrared dynamics in these two models are equivalent and, as \eB\ — > oo, 
the models reduce to a continuum set of (1 + 1) -dimensional GN models. 

In this paper, as in Refs. @-f§, only the case of a homogeneous magnetic field has been considered. To extend 
the present results to inhomogeneous field configurations, we note first that the number of colors in the GN model is 
Y c = (tt/2C)N c , where the factor n/2C — a 2 \eB\/2-K is proportional to the local magnetic flux attached to each point 
in X2a;3~plane. Let us consider a magnetic field B(xjJ, directed in +xi direction but depending on the transverse 
coordinates. It is tempting to speculate that in this case, as |e.B|(x_i_) — * oo, the NJL model will be reduced to a set of 
the GN models with the number of colors N c (x±) = (ir /2C(x±))N c (where C(xj_) = A 2 /|ei3|(x^)), which is different 
in different points of the a^xa-plane. It would be interesting to check these speculations by studying the NJL model 
in inhomogeneous field configurations. 

The results of the present paper are in agreement with the general conclusion of Refs. that the catalysis 

of chiral symmetry breaking by a magnetic field is a universal, model independent effect. This catalysis may have 
possible applications to cosmology, particle physics and condensed matter physics 

A specific application of interest would be a linkage between catalysis of chiral symmetry breaking and the existence 
of very strong primordial magnetic fields in the early universe Jul . The results of the present paper may be especially 
relevant for cosmological models based on supersymmetric dynamics |j) . 

Also, the effect of the dimensional reduction by external fields may be quite general and relevant for multi- 
dimensional field theories. 
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APPENDIX: EFFECTIVE ACTION IN THE SUSY NJL MODEL 



In this Appendix, the effective action for the SUSY NJL model in a magnetic field is derived. 
Besides the expression for the fermion propagator in a magnetic field, which was used in Refs. we also need 

to know the expression for the propagator in a magnetic field for a charged scalar with the mass m — p Wn: 



D(x. y) = exp 
where the Fourier transform of D(x) is 
D(k) = 



-( x - y yA<f{x + y) 



D{x - y) , 



ds 



cosheBs 



exp 



s(p2 _ fe 2 + ki tan^ £ ) +fei 



(Al) 



(A2) 



The expression ( ffl| ) for the effective potential in the SUSY NJL model can be obtained in the same way as for the 
ordinary NJL model in Ref. M. Let us show that this expression is equivalent to expression (H). 
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Eq.([ll|) can be rewritten as 



V(p) 



Nr. 



2G 16tt 2 Z 4 



In 



7 + 2/((p/) 2 ) 







where 



Let us show that 



I((3) = / dse 



-0s 



cosh s — 1 1 
s 2 sinh s 2s 



13 > 0. 



1 



1 



(/8+l)/2 



I(/3) = 1 - In £ + - In/3 - - IiiStt 2 + 4 



(ia; In T(x). 



0/1 



Using the integral representation for generalized zeta function , we find 

oo 

I(/3) = lim I(p,n) = lim / dss^e~ Ps * ' ' 



= lim T^) 



>• sinh i 



In order to find this limit, we need the following identity: 

92 - qi 



d_ 

dz 



C(z,q2) - ({z,qi) 
To derive it, we use the relation JlOj 



(<fe + qi - 1 - ln27r) + / dglnr(g) 



9g 



Differentiating it with respect to z and integrating over we obtain 

92 '/J 

CO*, 22) - C(z,qi) 



d_ 

dz 



J dqC{z + l,q) - z J dq^{z + l,q). 



Then, using the identities 



C(Q,q) = --q, 



we obtain Eq. (|A7|) . 

^From Eqs.(|A6D, (A7) and relations [p 



z=0 



lnr(g)--ln27r, 



iy + 2) 

(m + i)m' 



C(-i,«) = - 



q 2 - q 



6/ ' 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



we obtain equation (A5). 

Eqs.(|A3|) and (A5) lead to relation (|4j) for the effective potential. 

Let us show that, at weak coupling, the kinetic term in the effective action of the SUSY NJL model coincides with 
that in the ordinary NJL model. For this purpose, we shall prove that the contribution of the scalar fields in the 
kinetic term is suppressed. 



12 



By using the approach of Ref. j| , one finds that the functions f^ v and 



ft?" in the kinetic term 



to: 



Jl — \Jl If ex 
J 2. \J 2 ! j er 



2ip 2 



d A k - d 2 D(k) 



(27T)' 



-D(k) 



5|) are now equal 

(A12) 
(A13) 



where (f^ 1 ) j er are as given in 



(|l6|). The crucial point for us is that, while the fermionic contribution ( |28| ) in f? u 
is divergent as pi — > 0, the contribution of the scalars is finite in this limit. This conclusion, following directly from 
Eqs.(|A2|) and ( A13), reflects the fact that there is a gap in the spectrum of massless charged scalars in a magnetic field 
(see Eg. ([sij) ) . Since, at weak coupling, the parameter pi is exponentially small, we conclude that the contribution of 
the scalars is indeed suppressed at small g. 
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FIGURE CAPTIONS 

Figure 1. The curve of ntdynl = mdyn/ Vf^I as the function of the inverse coupling constant 1/g in the SUSY NJL 
model. 
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Figure 1 Elias et al Physical Review D 
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